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Abstract. - Electrical conductivity measurements on a 2D packing of metallic beads have 
been performed to study internal rearrangements in weakly perturbed granular materials. Small 
thermal perturbations lead to large non gaussian conductance fluctuations. These fluctuations 
are found to be intermittent in time, with peaks gathered in bursts. The distributions of the 
waiting time At between two peaks is found to be a power-law At~ ( - 1+at ' 1 inside bursts. The 
exponent at is independent of the bead network, the intensity of perturbation and external 
stress. These bursts are interpreted as the signature of individual bead creep rather than 
collective vault reorganisations. We propose a simple model giving at = 1 — £, where ( is the 
roughness exponent of the surface of the beads. 



Introduction. - Granular materials present interesting and unusual properties M. For 
example, photoelastic visualizations of static confined granular packing under external stress 
have provided evidence of large inhomogeneities in the distribution of contact forces between 
grains 0,^], giving rise to strong force chains (vaults or arches) extending on a scale much 
larger than the size of an individual grain. Recent experiments have revealed the sensitivity 
of sound transmission in granular media to small perturbations (such as thermal expansion of 
the grains) [Q . This sensitivity has been interpreted using a simple theoretical model of vault 
formation, the 'Scalar Arching Model' |^|. This model predicts that large rearrangements 
of the force chains network can indeed occur, even for small external perturbations. This 
sensitivity tends to disappear when the media is subjected to a strong external stress ||. 

However, both photoelastic visualization and sound transmission techniques require large 
external stresses and can not give information on a weakly confined granular media where 
interesting signatures of fragility |?j can be expected. We propose here to study this fragility 
through conductivity measurements. Experiments on the electrical properties of granular me- 
dia have been performed in the past : Branly reported || the influence of an electromagnetic 
wave on the electrical resistance of a granular packing. Giraud et al. ||, Marion et al. JlO| 
and Gervois et al. pd| ] have performed measurements of the electrical conductivity in order to 
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reveal the influence of the contact network on the mechanical properties of the media. More 
recently, Vandewalle et al. |l2| reported power-law distributions of conductivity changes, and 
argued that this was a signature of large rearrangments of stress paths. 

The aim of this work is to investigate whether force chain rearrangements can be observed 
through the variation of the electrical conductance of a 2D bead packings under small thermal 
perturbations. Large intermittent fluctuations are indeed observed. They appear in bursts 
that can be characterized via a statistical analysis of the waiting time separating two fluctu- 
ation peaks. The microscopic origin of these electrical fluctuations is then investigated. We 
argue that most of the interesting statistical features are not related to large force chains rear- 
rangements but to individual microcontacts between two beads that rearrange. The observed 
non trivial statistics might be related to the self-affine roughness of the surface of the beads. 

Experiment. - The experimental set-up is illustrated on fig. |l|. It consists mainly of a 2D 
packing of stainless steel beads of diameter a = 6mm± 10 /im confined between two Plexiglas 
plates (thickness 2 cm) . The use of metallic beads of millimetric size allows one to have a 
good control of geometrical and mechanical properties. We can also change the number of 
grains in the sample to separate local behaviour from collective behaviour. Beads are carefully 
arranged to form triangular compact packing. This situation of maximal compacity can easily 
be reproduced before each experiment. Even if the piling can be viewed as regular, the 
weak polydispersity and the solid friction between grains is sufficient to entail disorder in the 
network of contact forces []l3| . The size of the packing is W x H where W is its width and H 
its height in bead size units. Thermal perturbations are induced by a 75 W lamp standing at 
a distance d — 10 cm. The lamp increases the temperature of the beads and of the Plexiglas 
plates by typically 3 °C. This leads to an expansion of both the beads and the plates: the 
thermal expansivity for beads is 2 10 _5 K _1 and 7 10 _5 K _1 for Plexiglas. The bead network 
can also be vibrated via a buzzer placed against the Plexiglas plate. The packing is connected 
to a 9 V battery. A resistor R\ in series insures that the current crossing a contact between 
any couple of beads is much smaller than 40 mA. Preliminary studies performed on 2 beads 
indeed show that, as long as this current is below 40 mA, micro- welding is prevented and 
the contacts behave reversibly like an ohmic resistance. Voltages Vo, Vi and V2 (see fig. |l|) 
are recorded via a 12 bytes AD converter: Vq is the emf of the battery, Vi is proportional 
to the current crossing the packing and V2 represents the fast fluctuating part (faster than 
t = RC < 0.02 s) of V\. The following quantities can then be deduced: 

• the conductance g p of the packing. For a 33 x 52 vibrated packing, g p ~ 0.02 f2 _1 . 

• the fluctuation of conductance Ag p defined as the variation of g p between 2 acquisition 
points. Typically, Ag p reaches values up to 5.10~ 3 SI -1 to be compared to the detection 
threshold discussed below, which is typically around 5.10 _6 f2 _1 . 

A series of experiments was performed with the bead packing replaced by a 50 resistor 
in order to determine the level of noise of the acquisition system. This is mainly dominated by 
digital errors on V2 equal to 5 nV. A threshold is then defined: Ag p is considered as relevant 
when the corresponding AV2 is superior to 50 jiW . Sampling rates / from 0.05 Hz up to 10 kHz 
have been tried: / = 500 Hz is sufficient to separate most successive events (cf fig. [2]). For each 
experiment, data is recorded during a period of 20 minutes. The procedure is the following : 
(i) the packing is vibrated, (ii) the acquisition is started, (iii) the light is turned on 1 minute 
after (turned off in experiments with packing prepared with light switched on). 

Results. A packing, initially vibrated, keeps a constant conductance as long as its 
temperature does not change. As soon as it is thermally perturbed (light turned on or off), g p 
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starts to decrease and to fluctuate. Let us focus firstly on the fluctuations Ag p for a typical 
experiment (see fig. ||) on a 33 x 52 packing. Fluctuation peaks appear in an intermittent 
way. The distribution of Ag p is to a good approximation symmetrical and can be described 
(at least for small enough Ag p |) by the following power-law: 

P(\Ag p \) ~ |A ffp |-( 1 +« 9 ) (1) 

For this experiment, a g ~ 1, in agreement with the results reported in jl^]. For large \Ag p \, 
P(| A<7 P |) departs from this power-law behaviour. One can also study, for the same experiment, 
the temporal structure of the apparition of these peaks. The expanded view fig. |b of Ag p (t) 
presents a burst structure, which can be characterised by the distribution of the waiting time 
At between two successive events (independently of their sign). Two distinct regimes can be 
observed in fig. ||d: 

• The short At part of the probability density P(At) (corresponding to events within the 
same burst) decays as a power law with At: 

P(At) ~ At~ (1+at ) (2) 

The exponent is found to be at — 0.6. 

• for larger At, P(At) again departs from this power-law behaviour. The cumulative 
distribution function F(At) is less noisy for large At and shown in the inset of fig. |^d. 
One clearly sees that F decays exponentially for large At, which is the signature of a 
Poisson flux. The bursts themselves therefore appear as completely independent events, 
whereas the events inside bursts are strongly non-Poissonian. 

In order to determine the origin of these electrical fluctuations, several experiments have 
been performed with different packing: (1) regular 33 x 52 packing, (2) disordered packing 
of 1690 beads, (3) regular 33 x 32 packing, (4) regular 33 x 16 packing, (5) regular 33 x 2 
packing. Finally, (6) is a regular 33 x 2 packing where all metallic beads have been replaced 
by insulating glass beads except 3 beads forming a triangle in the middle. For all these 
experiments, P(\Ag p \) decays as a power-law but the exponent a g depends on the bead 
network geometry, and lies between 0.8 and 1.6 (see fig. ||a). On the contrary, the power-law 
decay of P(At) does not depend on the packing, and leads to a t ~ 0.6 (see fig. ||d). For the 
packing of 3 beads, the power-law behaviour actually extends over the three decades of the 
time distribution: only few bursts are present. The exponent is found to be at = 0.6 ± 0.05. 
When the number of beads increases, the value of At above which P(At) starts to deviate 
from a power-law decreases. It can be interpreted by the increase of independent burst number 
making the inter-burst waiting time mix more and more with the intra-burst statistics. The 
intensity of the thermal perturbation has been changed by modifying the distance of the lamp 
(d — 5 cm, 10cm and 20cm) without affecting the value of at- Moreover, for packing (5), 
the weight of the upper electrode has been varied (M — 108 g, M = 206 g, M = 304 g and 
M = 404 g). The number of events decreases, but -P(At) remains a power-law with the same 
exponent at = 0.6. 

Discussion. The origin of these non trivial electric fluctuations should be found in 
the geometry of electrical paths. A natural idea would be to relate these to force chains 
rearrangements. Indeed the Scalar Arching Model (SAM) introduced by two of us || suggests 
the existence of large force chains rearrangements in a packing subjected to small thermal 
perturbations. The SAM predicts a broad distribution of the apparent weight fluctuation AW a 
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measured at the bottom of the packing. This broad distribution of AW a could in principle 
be related to the observed broad distribution of Ag p . Moreover, the model also predicts 
avalanching contact reorganisations whose electrical signature could indeed be the bursts of 
fluctuations peaks seen in the thermally perturbed packing, each Ag p peak corresponding 
to the creation or the breakdown of a contact somewhere in the packing. However, in this 
scenario, the intra-burst structure should become less and less pronounced when the number 
of beads decreases, and should disappear completely for a packing of 3 beads. This is clearly 
not the case (see fig. ||b). 

The fact that the power-law nature of the distributions is independent of the geometry 
of the packing, and still holds for three beads, suggests that the origin of these electrical 
fluctuations is local. The direct imaging of the surface of the beads with an Atomic Force 
Microscope allowed us to measure the average roughness R of the beads. We have found 
R = 90 nm, to be compared to Hertzian deformation 5 in the packing. For contact forces on 
the order of 10~ 2 Newton, 6 ~ 10 nm which is smaller than R. The electrical contact between 
2 beads is thus completely dominated by the surface roughness. We propose then the following 
scenario to explain the observed behaviour of g p and Ag p in our experiment: when the packing 
is vibrated, the surface roughness of the two surfaces in contact 'adapt' to one another and 
the effective contact surface is high. Therefore, g p is maximal. As soon as the packing is 
heated, beads expand on the order of 0.5 /im. This thermal expansion is hindered by the 
roughness induced friction between beads. Therefore, most of the time the contact between 
two beads does not change, until the accumulated stress due to the thermal expansion is 
sufficient to 'unpin' the rough surfaces and make them slip over each other. This leads to the 
burst structure of the signal. Inside each burst, individual peaks of Ag p correspond to the 
creation or the loss of a micro-contact between the two slipping beads. These micro-contact 
avalanches globally lead to a smaller effective contact surface, thereby making g p smaller. 

Simple models of pinning [|l4|,[l5| can lead to a power-law distribution of the 'trapping' 
time between successive events; however, in these models the exponent at is found to depend 
on the external driving force. The fact that at is independent of both the strength of the 
perturbation and of the external stress suggests a purely geometrical interpretation: if the bead 
surface can be considered to be self-affine with a Hurst exponent the distance A£ separating 
two successive 'spikes' in a given direction is power- law distributed with an exponent a# = 1— C 
on the interval [01,02] where a± is a microscopic length below which the surface is flat, and 
<Z2 is at most the diameter of an Hertzian contact: 02 — 1 ^m. The typical speed V with 
which a bead thermally expands is 1 /zm per second. If this speed is constant in time, the time 
between successive micro-contact closing/opening is therefore a power law distribution with 
an exponent at — 1 — C m the interval [a\/V, (12/V — 1 s], which is precisely the experimental 
interval of time scales. We have determined ( from AFM frames via the moving average 
method the moving average of a series z{x) on a length interval L is defined as : 



For a self affine signal with an Hurst exponent of Q, the density p of crossing point between 
two moving averages on length scales L\ and L2 > L\ is found to scale as: 




L-l 



(3) 



i=0 



p~—[(AL)(l-AL)]S-i 



(4) 



where AL — (L2 — L\)jL<2, 
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AFM frames are stored as 256 x 256 matrices. Li is fixed to 128 and p(AL) has been 
determined for rows and columns averaging (see fig. [I]). In both case, we find £ = 0.35 ± 0.05, 
and therefore 1 — ( ~ 0.65. This value is consistent with at = 0.6. 

Conclusion. - In summary, a series of electrical measurements has been performed to 
test the SAM predictions of large force chains rearrangements in a 2D metallic beads packing 
under small perturbations. As soon as the packing is heated or cooled, its conductance starts 
to decrease and to fluctuate. These fluctuations are intermittent and gather in bursts as 
predicted by the SAM. However, the waiting time distribution inside bursts are independent 
of beads network, perturbation intensity and external applied stress. Consequently, the origin 
of these fluctuations should be found in local microcontact rearrangements at each beads 
rather than collective vaults rearrangements, as suggested in Jllfl . However, some information 
on these collective rearrangments might lie in the deviation from a pure power-law distribution 
which is clearly observed on figure 4 when one increases the number of grains. These collective 
effects might also be needed to account for the variation of the exponent a g with the geometry. 
Work in this direction is underway. 
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Fig. 1 - experimental set-up 
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Fig. 2 - Analysis of Ag p for a 33 x 52 regular packing heated with the lamp, (a) Temporal evolution 
of the g p . (b) Temporal evolution of Ag p . Inset, expanded view on a few seconds scale (c) Log-log 
plot of the distribution P(\Ag p \). The straight line is a power law fit using Eq. [l] with a g = l.(d) 
Log-log plot of the probability density function P(At). The straight line is a power law fit using 
Eq. |i| with OL t = 0.61. The inset shows the cumulative distribution which decays asymptotically as 
an exponential. 
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Fig. 3 - Log-log plot of (a) P(\Ag p \ and (b) P(At) for packing (1), (2), (3), (4), (5), and (6) from 
bottom to top. Data have been shifted for clarity. The straight lines are power-law fits. The exponent 
at has been fixed for all P(At) to at = 0.6 minimising errors on the fit for the three bead packing 
(6). a g is equal, from bottom to top, to 1.0, 1.1, 0.8, 1.6, 0.9, and 0.8 respectively. 
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Fig. 4 - (a) AFM picture of the bead surface (b) density p as a function of the relative difference 
AL with L2 = 128. o corresponds to row averaging and + correspond to column averaging. Inset : 
log 10 (p) as a function of log 10 (AL(l — AL)). The straight line corresponds to a fit using Eq. 



